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We show that the appHcation of a novel gauge invariant truncation scheme to the Schwinger- 
Dyson equations of QCD leads, in the Landau gauge, to an infrared finite gluon propagator and a 



4^ divergent ghost propagator, in qualitative agreement with recent lattice data. 
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Introduction - The infrared sector of Quantum Chromodynamics (QCD) [1] remains 
largely unexplored, mainly due to the fact that, unlike the electroweak sector of the Standard 
Model, it does not yield to a perturbative treatment. The basic building blocks of QCD are 
the Green's (correlation) functions of the fundamental physical degrees of freedom, gluons 
and quarks, and of the unphysical ghosts. Even though it is well-known that these quantities 
are not physical, since they depend on the gauge-fixing scheme and parameters used to 
quantize the theory, it is widely believed that reliable information on their non-perturbative 
structure is essential for unraveling the infrared dynamics of QCD [2] . 

The two basic non-perturbative tools for accomplishing this task are (z) the lattice, where 
space-time is discrctizcd and the quantities of interest are evaluated numerically, and {ii) 
the infinite set of coupled non-linear integral equations governing the dynamics of the QCD 
Green's functions, known as Schwinger-Dyson equations (SDE) [3-5]. Even though these 
equations are derived by an expansion about the free-field vacuum, they finally make no 
reference to it, or to perturbation theory, and can be used to address problems related to 
chiral symmetry breaking, dynamical mass generation, formation of bound states, and other 
non-perturbative effects [1]. While the lattice calculations are limited by the lattice size used 
and the corresponding extrapolation of the numerical results to the continuous limit, the 
fundamental conceptual difficulty in treating the SDE resides in the need for a self-consistent 
truncation scheme, i.e., one that does not compromise crucial properties of the quantities 
studied. 

It it generally accepted by now that the lattice yields in the Landau gauge (LG) an 
infrared finite gluon propagator and an infrared divergent ghost propagator. This rather 
characteristic behavior has been firmly established recently using large-volume lattices, for 
pure Yang-Mills (no quarks included), for both SU{2) [6] and SU{?>) [7]. To be sure, lat- 
tice simulations of gauge-dependent quantities are known to suffer from the problem of the 
Gribov copies, especially in the infrared regime, but it is generally believed that the effects 
are quantitative rather than qualitative. The effects of the Gribov ambiguity on the ghost 
propagator become more pronounced in the infrared, while their impact on the gluon prop- 
agator usually stay within the statistical error of the simulation [8] . In what follows we will 
assume that in the lattice results we use the Gribov problem is under control. 

In this article we show that the SDEs obtained within a new gauge-invariant trunca- 
tion scheme furnish results (in the LG) which are in qualitative agreement with the lattice 
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data. As has been first explained in [9], obtaining an infrared finite result for the gluon 
self-energy from SDEs, without violating the underlying local gauge symmetry, is far from 
trivial, and hinges crucially on one's ability to devise a self-consistent truncation scheme 
that would select a tractable and, at the same time, physically meaningful subset of these 
equations. To accomplish this, in the present work we will employ the new gauge-invariant 
truncation scheme derived in [10], which is based on the pinch technique [9, 11] and its 
correspondence [12] with the background field method (BFM) [13]. 

SDEs in the gauge-invariant truncation scheme - The gluon propagator A^j,(gf) in the 
covariant gauges assumes the form 



(1) 



where ^ denotes the gauge-fixing parameter, P^uiQ) — 9ni> — Q^Qu/q^ is the usual transverse 
projector, and, finally, A~^{q^) — + iU{q^), with n^,,(?) = P f^i,{q) U{q'^) the gluon self- 
energy. In addition, the full ghost propagator -D(p^) and its self-energy -L(p^) are related by 
iD~^{p'^) = — iL{p'^). In the case of pure (quarkless) QCD, the new SD series [10] for the 
gluon and ghost propagators reads (see also Fig. 1) 

[1 + G'(g2)]2 

iD-^p-") ^p' + tX f T^A,,{k)T''{p, k)D{p + k) , 
.Ik 

iA^.(g) - A / H^^^^D{k + g)A^<^(A;) H,.{k, q) , (2) 
Jk 

where A = g'^Cx, with Ca the Casimir eigenvalue of the adjoint representation [Ca = N 
for SU{N)], and = //^^(27r)~'^ / d'^k, with d — A — e the dimension of space-time. is 
the standard (asymmetric) gluon-ghost vertex at tree-level, and the fully-dressed one. 
G{q^) is the g/j^i, component of the auxiliary two-point function A^,^(g'), and the function 
Hcru is defined diagrammatically in Fig. 1. H^i, is in fact a familiar object [1]: it appears 
in the all-order Slavnov- Taylor identity (STl) satisfied by the standard thrcc-gluon vertex, 
and is related to the full gluon-ghost vertex by q'^H(^^{p,r,q) = —iT,^{p,r,q); at tree-level. 

When evaluating the diagrams (cj) one should use the BFM Feynman rules [13]; notice 
in particular that {%) the bare three- and four-gluon vertices depend explicitly on 1/^, {ii) 
the coupling of the ghost to a background gluon is symmetric in the ghost momenta, {iii) 




FIG. 1: The new SDE for the gluon-ghost system. Wavy hnes with white blobs are fuh gluon 
propagators, dashed hnes with white blobs are full- ghost propagators, black blobs are full vertices, 
and the grey blob denotes the scattering kernel. The circles attached to the external gluons denote 
that, from the point of view of Feynman rules, they are treated as background fields. 

there is a four-field coupling between two background gluons and two ghosts. Thus, for the 
gluonic contributions we find 

= ^ / f ^«/3A"^(A;)f ,,,A'^'^(fc + q) , 
^ Jk 

{a2),u = -i\g,. jA>;{k) - - l) j\u{k), (3) 

with r'^°'%q,pi,p2) = T'^°'%q,pi,p2) + {p2g'^'^—Pig'^^)^^\ ^f_iai3 the standard QCD three-gluon 
vertex, and F^q,^ is the fully-dressed version of F^^^. For the ghost contributions, we have 
instead 

(as)^, = -\ [ f^D{k)D{k + q)f, , 
Jk 

{ai)t,u = '2i\gf,u / D{k) , (4) 
Jk 

with T^{q,pi,p2) = {p2 —pi)ij., and its fully-dressed counterpart. Due to the Abelian 
all-order Ward Identities (Wis) that these two full vertices satisfy (for all namely 
l^^fia/s = "^^afii^ + l) ~ ^A^p{k) q^T ^ = iD~^{k + q) — iD~^{k), one can demonstrate 
that g'^[(ai) + (02)]^,^ = and q^Ka^) + (04)]^^ = [14]. 

For the rest of the article we will study the system of coupled SDKs (2) in the LG (^ = 0), 
in order to make contact with the recent lattice results of [6, 7]. This is a subtle exercise, 
because one cannot set directly ^ = in the integrals on the rhs of (3), due to the terms 
proportional to 1/^. Instead, one has to use the expressions for general ^, carry out explicitly 
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the set of cancellations produced when the terms proportional to ^ generated by the identity 
k^An^{k) — —i^kiy/k"^ are used to cancel 1/^ terms, and set ^ = only at the very end. It is 
relatively easy to establish that only the bare part T^a0 of the full vertex contains terms that 

diverge as ^ — ^ 0. Writing T,yaf3 = ^uai3 + ^iyai3, wc thus have that Kj,„^ is regular in that 
limit, and we will denote by K.^ai3 its value at ^ = 0. Introducing A^^j^{q) = P^j,(g)A(g^), we 
get 

where L/,Q^ = r^a/3+K^a/3 satisfies the WI = Pap{k + q)A~^{k + q) - Pap{k)A~^{k). 

Contracting the Ihs of (5) by q^ one can then verify that it vanishes, as announced. 

Next, following standard techniques, we express L^q^ and as a function of the gluon 
and ghost self-energy, respectively, in such a way as to automatically satisfy the correspond- 
ing Wis. Of course, this method leaves the transverse {i.e., identically conserved) part of 
the vertex undetermined. The Ansatz we will use is 



^ml3 = + [^apik + q) - ^a^ik)] , 

= f,-t^^[L{k + q)-L{k)] , (6) 

whose essential feature is the presence of massless pole terms, l/g^. Longitudinally coupled 
bound-state poles are known to be instrumental for obtaining A~^(0) ^ [15]; on the other 
hand, due to current conservation, they do not contribute to the ^'-matrix. For the conven- 
tional ghost-gluon vertex T^, appearing in the second SDE of (2) we will use its tree-level 
expression, i.e., T^, ^ — —p^. Note that, unlike T^, the conventional F,^ satisfies a STI 
of rather limited usefulness; the ability to employ such a different treatment for F^, and 

without compromising gauge-invariancc is indicative of the vcrsatiUty of the new SD 
formalism used here. Finally, for Ha-,y wc use its trcc-lcvcl vahic, Hi^J. 

With these approximations, the last two equations of (2), together with (4) and (5), give 
(in Euclidean space) 



[l + G{q')rA-\q') ^q'-^ 
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Since [(ai) + (a2)]ftj/ and [(0,3) + (04)]^!^ arc transverse, in arriving at (7) wc have used [(0-1) + 
(02)]/^!/ = Tr[(ai) + {a2)]Piuu{q) and [(03) + (04)]^!, = Tr[(a3) + {a4)]Pf,„{q), substituted into 
(2), and then equated the scalar co-factors of both sides. Thus, the transversahty of the 
answer cannot be possibly compromised by the ensuing numerical treatment (e.g. hard 
ultraviolet cutoffs), which may only affect the value of the co-factor. 

Numerical results - Before solving numerically the above system of integral equations, 
one must introduce renormalization constants to make them finite. The values of these 
constants will be fixed by the conditions A~^(//^) — jj?, D~^[iJ?) — jj?, and = 0, 

with the renormahzation point jj? of the order of M|. It is relatively straightforward 
to verify that the perturbative expansion of (7) and (8) furnishes the correct one-loop 
results. Specifically, keeping only leading logs, we have l-\-G{q^) — 1 -\- ^^g"" ln(g^///^), 

ii|^ln(gV/x2)], where 



while D-\p'')=p^[l + ^^\n{p'/^l'')] and /^-\q^) = q^[l ^ 
as = g^/^n. 

The crux of the matter, however, is the behavior of (7) as q^ 
of the gluon propagator is observed. In this limit, Eq.(7) yields 

[1 + ^(0)12' 

15 



/ Aik) - I I eA\k), 
Jk ^ Jk 

I D{k)+ I k^D'^ik). 

Jk Jk 



0, where the "freezing" 

(9) 
(10) 

(11) 



Perturbatively the rhs of Eq.(9) vanishes by virtue of the dimensional regularization result 
/,^ = 0n = 0,1,2,... which ensures the masslessness of the gluon to all orders. How- 
ever, non-perturbatively A~^(0) does not have to vanish, provided that the quadratically 
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divergent integrals defining it can be properly regulated and made finite, without introducing 
counterterms of the form mQ(A^y)A^, which are forbidden by the local gauge invariance of 
the fundamental QCD Lagrangian. It turns out that this is indeed possible: the divergent 
integrals can be regulated by subtracting appropriate combinations of "dimensional regular- 
ization zeros". Specifically, as we have verified explicitly and as can be clearly seen in Fig.2 
(left panel), for large enough /c^ the A(A;^) goes over to its perturbative expression, to be 
denoted by Apert(A;^); it has the form Apert(A;^) = '^n=o ^n^^^^, where the coefficient a„ are 
known from the perturbative expansion. For the case at hand, measuring in GeV^, using 
fj, 100 GeV and as{fJ.) = 0.1, after inverting and re-expanding the A~^(/c^) given below 
Eq.(8), wc find oq ~ 1.7, ai ^ —0.1, 02 ~ 2.5 x 10~^. Then, subtracting Apert(/^^) = 
from both sides of Eq.(lO), we obtain the regularized T^^^ given by (/c^ = y) 



A similar procedure can be followed for Tc (sec below). The obvious ambiguity of the 
regularization described above is the choice of the point s, past which the two curves, A(y) 
and Apert(y), are assumed to coincide. 

Ideally, one should then: (i) solve the system of integral equations under the boundary 
condition A(0) = C, where C is an arbitrary positive parameter; (ii) substitute the solu- 
tions for A(g) and D{q) in the (regularized) integrals on the rhs of (9), together with the 
obtained value for G'(O), and denote the result by Aj:pg(0); (m) check that the self-consistency 
requirement Aj:pg(0) = is satisfied; if not, (iv) a new C must be chosen and the proce- 
dure repeated. In practice, due to the aforementioned ambiguity, we cannot pin down A(0) 
completely, and we will restrict ourselves to providing a reasonable range for its value. 

We have solved the system for a variety of initial values for C, ranging between 
1 — 50GeV^^, and obtained from (12) the corresponding A~g(0). On physical grounds 
one docs not expect the perturbative expression Apcrt(A:^) to hold below 5 — lOGeV^, and 
therefore, when computing A~g(0), s should be chosen around that value. For values of C 
between 10 — 25GeV~^ the corresponding A~g(0) can be made equal to by choosing 
values for s within that (physically reasonable) range. For example, for C = 14.7GeV-^ the 
value of the lattice data at the origin, we must choose s 10 GeV^. The solutions for A.{q), 
D{p), and 1 -|- G{q) obtained for that special choice, C — 14.7 GeV~^, are shown in Fig.2. 
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FIG. 2: Left panel: The gluon propagator obtained from the solution of the SDE system (blue 
continuous line) compared to the lattice data of [7] ; the red dashed line represents the perturbative 
behavior. In the inset we show the function 1 + G{q^) (blue continuous line) and its perturbative 
behavior (red dashed line). Right panel: The ghost propagator obtained from the SDE system 
(blue continuous line), the one-loop perturbative result (red dashed line), and the corresponding 
lattice data of [7]. In the inset we show the function p'^D{p^) from the SDE. 

In order to enforce the equality A~g(0) = C^^ for higher values of C one must assume the 
validity of perturbation theory uncomfortably deep into the infrared region; for example, for 
C = 50 GeV~^ one must choose s below IGeV^. We emphasize that the non-perturbative 
transverse gluon propagator, being finite in the IR, is automatically less singular than a 
simple pole, thus satisfying the corresponding Kugo-Ojima (KO) confinement criterion [16], 
essential for ensuring an unbroken color charge in QCD [17]. Note that for < 10 GeV^ 
both gluon propagators (lattice and SDE) shown in Fig.2 may be fitted very accurately 
using a unique functional form, given by A^^(g^) = a + b {q'^Y^^ . Specifically, [measuring 
in GeV^ and the per degrees of freedom], the lattice data are fitted by a = 0.07, b = 0.15, 
and c = 2.54 (x^ ~ 10~^), while our SDE solution is described setting a = 0.07, b = 0.77, 
and c = 2.01 (x^ ~ 10"^). 

Let us now consider the ghosts. The D{p'^) obtained from the ghost SDE diverges at 
the origin, in qualitative agreement with the lattice data. From the SDE point of view, 
this divergent behavior is due to the fact that we are working in the LG and the vertex 
Fj, employed contains no poles, as suggested by previous lattice studies [18]. The 
rate of divergence of our solution is particularly interesting, because it is related to the 
KO confinement criterion for the ghost [16], according to which the non-perturbative ghost 
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propagator (in the LG) should be more singular in the infrared than a simple pole. Motivated 
by this, we proceed to fit the function p^D{p^) [see inset in right panel of Fig.2]. First we use 
a fitting function of the form p^Dip^) = Ci{p'^)~^ {p^ in GeV^); a positive 7 would indicate 
that the SDE solution satisfies the KO criterion. Our best fit, valid for < 10, gives the 
values 7 = 0.02 and Ci = 1.30, which lead to a ~ 10"^. Interestingly enough, an even 
better fit may be obtained using a qualitatively different, physically motivated functional 
form, namely p^D(p^) — Ki — K2\n{p^ + K3) (with % acting as a gluon "mass"). Our best 
fit, valid for the same range, gives k-i — 1.3, K2 — 0.05, and = 0.05, with ~ 10~^. This 
second fit suggests that p^D{p^) reaches a finite (positive) value as p^ 0. Even though 
not conclusive, our fitting analysis seems to favor a ghost propagator displaying no power- 
law enhancement, in agreement with recent results presented in [19]; clearly, this question 
deserves further study. 

Turning to the tadpole contributions Tc of (11), the subtraction of = Xc regularizes 



Tc, yielding a rather suppressed finite value for TJ^s por example, using the first ghost fit, 
we get (s' 1 GeV^) 



which is numerically negligible. 

Discussion - The present work has focused on the derivation of an infrared finite gluon 
propagator from a gauge-invariant set of SDEs for pure QCD in the LG, and its comparison 
with recent lattice data. Following the classic works of [15], the finiteness of the gluon prop- 
agator is obtained by introducing massless poles in the corresponding three-gluon vertex. 
The actual value of A~g(0) has been treated as a free parameter, and was chosen to coincide 
with the lattice point at the origin. The curves shown in Fig.2 were then obtained dynami- 
cally, from the solution of the SDE system, for the entire range of momenta. Comparing the 
solution for the gluon propagator with the lattice data we see that, whereas their asymptotic 
behavior coincides (perturbative limits), there is a discrepancy of about a factor of 2-2.5 in 
the intermediate region of momenta, especially around the fundamental QCD mass-scale 
[refiected also in the different values of the two sets of fitting parameters (a,6, c)]. In the 
case of the ghost propagator the relative difference increases as one approaches the deep 
infrared, given that both curves diverge at a different rate. These discrepancies may be ac- 
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counted for by extending the gluon SDE to include the "two-loop dressed" graphs, omitted 
(gauge- invariantly) from the present analysis, and/or by supplying the relevant transverse 
parts of the vertex given in (6). We hope to be able to make progress in this direction in 

the near future. 

In our opinion, the analysis presented here, in conjunction with the recent lattice data, 
fully corroborates Cornwall's early description of QCD in terms of a dynamically generated, 
momentum-dependent gluon mass [9] . In this picture the low-energy effective theory of QCD 
is a non-linear sigma model, known as massive gauge-invariant Yang-Mills, obtained from 
the generahzation of Stiickelberg's construction to non-Abehan theories [20]. This model 
admits vortex solutions, with a long-range pure gauge term in their potentials, which endows 
them with a topological quantum number corresponding to the center of the gauge group 
[Zisf for SU{N)], and is, in turn, responsible for quark confinement and gluon screening 
[21, 22]. Specifically, center vortices of thickness ~ m~^, where m is the induced mass of the 
gluon, form a condensate because their entropy (per unit size) is larger than their action. 
This condensation furnishes an area law to the fundamental representation Wilson loop, thus 
confining quarks. On the other hand, the adjoint potential shows a roughly linear regime 
followed by string breaking when the potential energy is about 2m, corresponding to gluon 
screening [23]. 
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